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Abstract 

A relativistic and non-local model for the coherent photoproduction of pions 
on spin-zero nuclei is presented. The production operator is derived from 
an effective Lagrangian model that contains all the degrees of freedom known 
to contribute to the elementary process. Using the framework of DWIA, the 
matrix elements of this operator are evaluated using relativistic bound state 
wave functions. Final state interactions are accounted for via a pion- nucleus 
optical potential. The effects of in-medium modifications of the production 
operator are discussed. We give results for 12 C and 40 Ca, and compare our 
calculations to the data of the A2 collaboration at MAMI. 

I. INTRODUCTION 

Photonuclear reactions offer a unique possibility to test our understanding of hadrons in 
vacuum as well as within a nuclear environment. For elementary reactions like the photopro- 
duction of pions on nucleons, relativistic models based on an effective Lagrangian approach 
are well established by now (see for example (ij and references therein). Starting from the 
elementary process, numerous studies treat the photoproduction of pions on nuclei. Among 
the large number of possible reactions, the coherent photoproduction plays a special role. 
Since in this case the nucleus remains in the ground state after the reaction process, the 
amount of further theoretical assumptions is smaller than e.g. in the case of reactions like 
(7, ttN) or completely inclusive processes like photoabsorption on nuclei. In this sense the 
coherent photoproduction of pions is the cleanest test for any model that treats photopro- 
duction of pions on nuclei. 
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Data on the coherent photoproduction of pions on nuclei were published in ||. The 
extraction of the coherent events in most of these experiments, however, was uncertain, so 
that comparison to these data is not very conclusive. More recently, the differential cross 
section was given for a photon energy of E 1 « 175 MeV in ||. It is, however, a special 
feature of the coherent process, that the differential cross section is mainly sensitive to the 



nuclear form factor (see Sec. |II Q ). Since the nuclear wave functions, and therefore the 



bound state properties of the nucleus, are taken as an input, all models yield more or less 
the correct shape of the differential cross section as can be seen in |3j . For a given energy they 
differ mainly with respect to the absolute value of the cross section. To further distinguish 
between different models, energy dependent data are needed. Thus the new data of the 
A2-collaboration at MAMI Q, which give the differential cross section at a fixed angle for 
a large energy range between 140 MeV and 430 MeV, are highly welcome. They enable us 
to test our calculations not only close to threshold, but over the entire A region. 

Part of the existing calculations for the coherent photoproduction of pions work in the 
delta-hole model B-4T3"f. These calculations treat the A and pion dynamics in a very sophis- 



ticated manner, however, non-resonant contributions are normally neglected (See, however, 



P Jll]Jl2]| ). While many other works use a local approximation, in [12| the non- locality of 
the process is treated explicitly. Another group of models is mainly used for the thresh- 
old region and employs the distorted wave impulse approximation (DWIA) fT4]-18|. These 



studies take into account the non-resonant contributions, but resort for technical reasons to 
non-relativistic and local approximations. 

In this paper we present a relativistic, non-local approach to the coherent photopro- 
duction of neutral pions on spin-zero nuclei in the framework of DWIA. The production 
operator is taken from an effective field theory, which describes the cross section for the 
production of neutral pions on nucleons well. The wave function of the outgoing pion is 
calculated using an empirical 7r-nucleus optical potential including s- and p-wave contribu- 
tions. The bound state wave functions are extracted from a relativistic equation of motion 
with empirically determined scalar and vector potentials. We thus include resonant and 
non-resonant contributions to the production operator in a consistently relativistic man- 
ner, without non-relativistic or local approximations. We also study the effect of medium 
modifications of the production operator, especially in the A contribution. This approach 
enables us to test our understanding of the relativistic structure of the nuclear ground state 
and the production operator, though we resort to empirical models for the nucleus and the 
pion-nucleus interaction. 

In the following section we first describe the elementary production operator we use. 
After some general remarks on the coherent photoproduction of pions, we then give the 
details of our model. Results for PWIA, as well as DWIA calculations are shown in Sec. [II. 
We finally discuss the effects of in-medium modifications of the production operator. In 
Sec. [TV] the results of our study are summarized. 



II. THE MODEL 
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A. The elementary operator 



The diagrams contributing to the elementary process of photoproduction of neutral pions 
on a single nucleon are shown in Fig. [1]: There are the direct and crossed nucleon graph, 
the direct and crossed A-graph and the graph containing the omega meson in the t-channel. 
There could in principle also be an analogous graph containing the rho meson, but due the 
smallness of the corresponding couplings it contributes only negligibly to the production of 
neutral pions JTJ. We use the following interaction terms: 

C iNN = _ e ^i(i + T3 ) 7jU ^A" 
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with F^ = d fJi A u — dvA^ and = d^uJv — d u u^, where A^uj^) denotes the photon (omega) 
field. 



For the A-propagator we take [19 



with 



G%(p) = i t + m ^ T , / s (2) 



A- = ( 9 --i^-^ pV -J-( 7V -^)) . (3) 



The width of the A is taken to depend on the energy: 

\Qo) \q 2 + c 



r W = r^(l)'(i±4)' (4) 
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with r o = 120 MeV and c = 0.5 GeV. q and q Q denote the momentum of the pion resulting 
from a decay of a A of a mass ^/s and of 1.232 GeV, respectively, in the rest-frame of the 
A. 

For the coupling constants we use the following values: 
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= 0.3028 


J 9ttNN 


= 0.97 




= 1.79 




= -1.91 
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(5) 



The coupling constants of the photon and the pion to the nucleon are well known. The 
value we use for g wNA corresponds to an on-shell decay width of the A of 120 MeV, so that 
it is consistent with T in Eq. (Q). g has been determined from the decay width of an 
omega meson into a pion and a photon [|TJ. For the couplings of the omega to the nucleon 
different values are used in the literature; the values given in (^]) are taken from [p0[| . The 
contribution of the omega graph to the photoproduction of pions or eta mesons is small, so 
that the uNN coupling constants cannot be determined by fitting the experimental data 
for these processes [^l],[22|]. Instead, they are taken from NN scattering |p3| |, or they are 



derived from the pNN parameters via SU(3) considerations [23|. The values for g v used 



in the literature for photoproduction processes range from 8 [24] to 17 fl25 
form factor is introduced at the uNN -vertex: 



F(t) 



A 2 - m 2 
A 2 -t 



In addition, a 



(6) 



with A 2 = 2 GeV 2 |21|. The sensitivity of our results to the omega parameters will be 
discussed in Sec. fU . 

Finally, g iNA and c in Eq. were determined by comparing the results of the present 
model to the data for the cross section of the photoproduction of neutral pions on protons 



As can be seen from Fig. 



the differential cross section can be well reproduced 
with the parameters given above. We did not use multipoles for the determination of these 
parameters, since multipoles are especially sensitive to unitarization effects, which are not 
included in the present model. 



B. Photoproduction on the nucleus 

After having specified the elementary production operator we now turn to the photo- 
production on the nucleus. The framework commonly used for exclusive processes like the 
coherent photoproduction is the DWIA. Here it is assumed that the production process in- 
volves only a single nucleon, while for the distortions of in- or outgoing particles interactions 
with the entire nucleus are taken into account. Thus many-body contributions to the pro- 
duction process are not taken into account. The single-particle production operator is taken 
from an elementary model and evaluated using relativistic bound state and scattering wave 
functions instead of plane waves. Since the nucleon participating in this reaction is bound 
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in a potential, it is off-shell, and the production operator must be evaluated for kinematical 
situations different from the one in the free case. 

For technical reasons a relativistic production operator is often simplified using on-shell 
relations |27]] or a non- relativistic reduction |25| , in order to make its evaluation easier. The 
resulting production operators are more or less equivalent on-shell, but the off-shell behavior 
is not the same as for the original operator. To avoid this problem, we take the production 
operator in its original form from the model described in the previous section, without 
rewriting it. Thus we use the natural off-shell dependence resulting from an effective field 
theory. 

A numerical complication arises from the fact that the production operator depends on 
the momentum of the struck nucleon, i.e. it is non-local. To circumvent complicated integra- 
tions, most DWIA calculations evaluate the production operator in a local approximation 
at some fixed effective nucleon momentum. For knock-out processes like (7, irN) the asymp- 
totic momenta of the outgoing particles can be used to estimate the momentum arguments 
of the production operator. In the case of coherent production, however, both the incoming 
and the outgoing nucleon are in a bound state, so that they do not have a well defined 
asymptotic momentum. The validity of the local approximation for the photoproduction 
of charged pions is is discussed in To avoid the uncertainties related to the local 

approximation, we evaluate the matrix elements of the production operator non-locally. We 
will discuss the local and the non-relativistic approximation further in Sec. |TTT| . 



C. The coherent process 



In the nuclear coherent photoproduction of pions, the nucleus remains in the ground 
state after the reaction. Since the pion in the final state has the quantum numbers 0~, the 



Lorentz invariant amplitude for the entire process must have the form [ 27] : 



T (A) 



(7) 



with 



T^ = e^K Pp q a A(s,t) . (8) 

k, p and q are the momenta of the incoming photon, the incoming nucleus and the outgoing 
pion, respectively. eW is the polarization vector of the photon and A(s, t) is a scalar function 
that contains the entire dynamics of the process and depends on the Mandelstam variables 
s and t. This can be shown to yield: 

E I T= ^ *L iL^Ocm I A(S, t) I 2 . (9) 

A 

k 2 cm and q 2 cm are the three momenta of the photon and the pion, 9 cm is the scattering angle 
and W is the total energy of the photon-nucleus system; all these quantities are taken in 
the cm-frame. Thus the well known sin 2 ^-dependence of the coherent photoproduction of 
pions results directly from the quantum numbers involved. The differential cross section in 
the cm-frame is then given by: 
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where M4 is the mass of the nucleus. 

Another important consequence of Eq. (|7|) is, that if one replaces the photon polarization 
by the photon momentum k^, vanishes. Thus the amplitude is gauge invariant 
from the very beginning, independent of the model used for the nuclear ground state or the 
production process. This is a special feature of the coherent photoproduction on spin-zero 
nuclei, which makes this reaction even more attractive from a theoretical point of view: In 
other reactions than the coherent one, like (7, ttN) or (e,e'N), the usual DWIA approach 
leads to a gauge- dependent amplitude f30H3~2|. Methods to restore gauge invariance lead to 
theoretical uncertainties |31| , which do not occur in the case of coherent photoproduction. 



We work in position space, since the bound state and scattering state wave functions can 
easily be obtained in a position space representation. In the approach described above the 
first diagram in Fig. [I] corresponds to the following non-local expression: 

4 A i=E /^w^)^*^^ ■ (ii) 



Here ip a is the wave function of the bound nucleon, (f)^ is the wave function of the photon, 
which can be assumed to be a plane wave and is the distorted wave function of the pion 
satisfying incoming boundary conditions p3| . T^nn and T^ NN are the vertices resulting 



from the coupling terms in Eq. (|1]) and G° N is the free nucleon propagator: 

. (12) 

The effect of including a dressed instead of a free nucleon propagator will be discussed in 



Sec. [Ill C[ E is naturally determined by energy conservation: 

E = E 1 + E a , (13) 

where E a is the total, relativistic energy of the bound nucleon. 

To evaluate (|ll]), we use partial- wave expansions for the wave functions as well as for 
the propagator [031 . Whenever derivatives appear in the vertices, we actually insert the 
derivatives of wave functions, so that we take into account the non-locality of the production 
operator. After inserting the partial-wave expansions into Eq. ( |Tl~D we can perform the 
angular integrations analytically, since they only involve spherical harmonics. The two 
radial integrations are evaluated numerically. 

The remaining graphs in Fig. |T] are treated analogously. For the A-propagator we treat 
the contraction of A MJ , from Eq. @ with the AA^7-vertex as one term and use the partial- 
wave expansion only for the remaining spin-1 part of the propagator: 



s - m\ + i-y/s T(s) 



by taking its Fourier transform analogously to Eq. (|T|). However, a closed expression of this 
Fourier transform, which is needed in order to use the formula for a partial wave expansion, 
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can only be obtained if the imaginary part of the denominator does not depend on the 
three-momentum p. Since s T(s) depends on p via s = p 2 Q — p 2 we evaluate this term at a 
value s Q (E 7 , E a ), which is taken to be the invariant mass of a system of the photon and the 
struck nucleon, averaged over the Fermi-sphere. Consequently, the A-width is a function of 
the energy of the incoming photon (i? 7 ) and the struck nucleon (E a ) . 

A further technical problem arises from the p^Pv term in Eq. (H), since its spatial compo- 
nents piPj contain second derivatives, an exact treatment of which would greatly complicate 
the calculations. We therefore approximate this term by replacing 

PiPj -> hkj , (15) 

where k is the three- momentum of the photon. Since in momentum space the three- 
momentum of the A-propagator is the sum of photon and nucleon momentum, this amounts 
to putting for this special term the three-momentum of the incoming nucleon to zero. Using 
this approximation, this term contributes only on the percent level, thus affecting our results 
insignificantly. 

The uiNN form factor in Eq. (|(|) is approximated by using the asymptotic momentum of 
the pion in t = (p 7 — p-w) 2 - In the relevant kinematic region, t is small anyway as compared 
to the cutoff A 2 . Consequently the dependence of F(t) on the pion distortions is weak and 
its main effect is a renormalization of the ujNN coupling. 

The angular integrations in Eq. ( p"T|) lead to rather complicated expressions for the ver- 
tices depending on the incoming and outgoing angular momenta. The sin^ dependence of 
the differential cross section is one test for their correct numerical implementation. As a 
further test we have checked that for T M from Eq. (§) our calculation yields: 

T o = 0;T i k i = , (16) 

which has to be fulfilled for each graph separately, as can be seen from the definition of this 
quantity. 

Even though we evaluate the amplitude in position space it is instructive to consider 
Eq. ( |ll~D in momentum space: 

T Nd ir= I TTTTz ^a(P + k- q) T nNN G N (E;p + k) T lNN tjj a (p) 

a occ.'' \^) 

^ p d^p -> -» 

= 12 / ^a(P + k - q) T Ndir (E;p, k,q) ip a (p) , (17) 

a occ. J 

where for the moment we assume the pion wave function to be a plane wave. If we now 
make the local approximation as described at the end of Sec. H B| by putting the momentum 



of the incoming nucleon in the argument of T^dir equal to a constant p Q , we can rewrite 
Eq. « 



T Ndi r ~ Tr (f Ndir (E; p Q , k, q) p A (k - g)) , (18) 
with the density matrix: 
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a occ. 



(2vr) ; 

d s xe ipx J- if) a (x) ^ a {x) = I d 3 x e ip_sr p A {x) . (19) 

a occ. 

The density matrix Pa{%) contains the complete information about the nuclear ground state. 
The scalar, vector and tensor densities are obtained by taking Tr(TpA) with T = l,7 ,cr *, 
respectively. This is in contrast to the approach used in |27| for the coherent photoproduction 



of //-mesons. There the production operator is rewritten using the free Dirac equation for 
in- and outgoing nucleons. As a consequence, the coherent photoproduction of ?7-mesons 
depends in p7| only on the tensor density of the nuclear ground state. 



One sees from Eqs. (18) and (0), that in the local approximation the differential cross 
section for coherent photoproduction at a fixed energy involves the Fourier transformation 
of a nuclear ground state density, i.e. a nuclear form factor. Although effects due to the 
non-locality of the production operator and the distortions of the pion wave function by the 
nucleus are superimposed, the differential cross section is dominated by the nuclear form 
factor. Hence the properties of the production operator itself can only be studied when the 
energy dependence of the cross section is considered. 



D. The nuclear wave function 

It is clear from these considerations, that it is crucial to describe the ground state prop- 
erties of the nucleus realistically. The wave functions in Eq. (|TT| ) are taken from a relativistic 
mean-field calculation using scalar and time-like vector potentials V s and V v , respectively: 



rn 



Vvlo - V s )i/j c 







(20) 



For the potentials V v and V s we assume a Woods-Saxon shape: 



V(r) 



1 + e" 



(r-TjA 1 / 3 ) 



v, s 



(21) 



The parameters we use for these potentials are given in Table |I[ They were determined such 
that the separation energies [Pjj, the root mean square radius of the charge density and the 
charge form factors of 12 C and 40 Ca [35| are well reproduced. The resulting charge form 
factors are shown in Figs. |3] and f|, together with the values extracted from experiment [35| . 



By using these wave functions, we neglect higher order correlations beyond mean field dy- 
namics. The question arises, to what extent the momentum structure of the nuclear ground 
state wave functions is adequately described by mean- field dynamics [j36,37||. The coherent 
production process leaves the nucleus in the ground state and probes, in the energy range 
considered here, primarily momenta inside the Fermi-sphere. Thus, a reliable description 
of nuclear form factors in a momentum range up to about twice the Fermi-momentum is 
most important. As seen in Figs. |3] and ^, the measured charge form factors are indeed well 
reproduced up to a momentum transfer of about 3 fm _1 , so that the relevant momentum 
range is obviously well described. Extended approaches (e.g. p8|,|39|) are in fact showing 
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that for spherical nuclei deviations from mean-field dynamics will become detectable only at 
rather high momenta corresponding to short-range processes which are not considered here. 
A mean field version of the Walecka model has recently been used to treat the coherent 
photoproduction of 77-mesons on nuclei [pT]]. 



E. The pion-nucleus interaction 



It is well known, that the coherent cross section depends strongly on the interaction of 
the produced pion with the nucleus. This is taken into account by using a distorted pion 
wave function in Eq. (|Tl"D, which is calculated from a position space optical potential. An 
optical potential as given in is used: 



2E w U opt (r) 



-47r[6(r) + B(r)} 
+4vrV{L(r)[c(r) + C(r)]}V 



-47T 



^V 2 c(r) + ^V 2 C(r) 



(22) 



where 



with 



b(r) 
c(r) 
B(r) 
C(r) 

L(r) 



Pi = l 



Pib p(r) 

Pi l Cop(r) 

p 2 B p 2 (r) 

1 + ^A[c(r) + C(r)f ~* 



m N ' 



P2 = l + 



2m 



N 



(23) 



E n is the total energy of the pion, is the mass of the nucleon and p is the nuclear density 
normalized to A, obtained by summing the single particle densities in the potentials Eq. fl2T|) 
over the occupied states. This optical potential contains s- and p-wave interactions of the 
pion via b(r) and c(r) and the so called true absorption via B(r) and C(r). The quantity 
L{r) comes in because of the Lorentz-Lorenz-Ericson-Ericson effect Jl9|]. The last term in 
Eq. ( |2"2"| ) results from the so called angle transformation jnj, which is a kinematical effect. 
The parameter A is a real constant, while the quantities b Q , c a , C a , and C Q are complex and 
energy dependent. In |^0| , ^l| these parameters are only given for pion kinetic energies up 
to 50 MeV. Since we want to compare our calculations to the A2 data, we need an optical 
potential for pion kinetic energies up to about 300 MeV. In order to have a consistent 
parameterization over this range of energies, we adopt the same form as in Eq. fl2"2"|) also for 
higher energies. The parameters were determined by fitting the elastic scattering data of 



pions on 12 C [12], including the nuclear Coulomb potential as done in |4(| . 

Doing so one is confronted with the fact that there is a large redundancy between some 
of the parameters, especially between b Q and B Q and between c Q and C Q [J43]. Consequently, 
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a naive fit of all parameters leads to unphysical values. We therefore put B = and C Q = 
and determine b Q and c a as complex, effective parameters. In Fig. |5| we show the cross 
sections for elastic pion scattering on 12 C in comparison to the experimental data. In order 
to have one more set of parameters at higher energies, we included elastic pion scattering 
on 16 at T w = 330 MeV into the fitting procedure. For T n = 100 MeV and T n = 157 MeV, 
we were not able to find a parameter set that led to agreement with the data beyond the 
second minimum of the differential cross section, without showing discrepancies at smaller 
angles. In these two cases we chose the optical potential parameters such that the data 
are reproduced well for angles up to the second minimum. In order to obtain the optical 
parameters as smooth functions of the pion energy, we have interpolated our fit results as 
shown in Fig. 

We expect that the empirical potential thus obtained describes the final state interactions 
of the pion in the case of 12 C sufficiently well in an energy range from threshold up to the A- 
resonance. The structure of this potential is different from potentials that were obtained from 
microscopic calculations [fi0| , fl4l1 , for which we also show results in Fig. [5| We will compare 
the effects of our potential in the coherent photoproduction to the effects of potentials that 
were extracted from microscopic models in Sec. [111 B| . 



III. RESULTS 
A. PWIA 



The model described in the previous section was used to calculate differential and total 
cross sections for the coherent photoproduction of pions on nuclei. In Fig. [7| and Fig. ^| we 
show our results for the plane wave impulse approximation (PWIA), i.e. with a free outgoing 
pion. As has been discussed at the end of Sec. [11 0, the total cross section shows the resonant 



character of the elementary production operator, while the shape of the differential cross 
section results from the nuclear form factor (multiplied with sin 2 #). In Fig. [7] we also 
show the contributions of the different diagrams. Clearly, the A-resonance dominates the 
process, but the other diagrams lead to sizable corrections that should not be neglected. 
Note the destructive interference between the Born terms and the other graphs. The omega 
contribution comes mainly from the vector coupling of the omega (Eq. ([!])), the tensor 
coupling contributes only negligibly. 

In order to compare to the non-relativistic, local approach used in the previous DWIA 
calculations, we performed a local calculation using the non-relativistic reduction of the 
direct A-graph as given in f25 |, which has also been used in fpT5 ,pT7 1 . The direct A-graph 



accounts for most of the total cross section (see Fig. [7]). In order to compare this to our 
result for the direct A-graph, we have used our set of parameters (Eq. ([5])) including the 
energy dependent A-width from Eq. ([|). The struck nucleon was assigned an energy equal 
to the average energy of a bound nucleon in 12 C (cf. Eq. (|T3"|)). The three-momentum ^ of 
the struck nucleon in the photon-nucleus cm-system was assumed to be 

k A 1 . — * . , . 

# = "I - "2^-(* - # , (24) 
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which is often used in local DWIA calculations (see e.g. fl7||). The nuclear form factor 
was taken to be the vector form factor resulting from our ground state calculation, which 
corresponds to the charge form factor shown in Fig. RL corrected for the electromagnetic 
form factor of proton and neutron. 

The result of this calculation for the direct A-graph is shown as the dotted curve Fig. |9| 
for 12 C together with the result of 'exact' calculations (solid curve). In 1)51,161 it is found 
that in a non-relativistic approach there are additional corrections due to the transformation 
between the photon-nucleus and the photon-nucleon cm-frame. In || it is shown that the 
inclusion of this transformation corresponds to multiplying the amplitude with a factor 
m[(m + E^y/l — v 2 ]" 1 , where m is the mass of the nucleon and v is the boost velocity 
from the photon-nucleus to the photon-nucleon cm-system. Including this factor in our 
non-relativistic, local calculation leads to the dashed curve in Fig. [|. It is an advantage of 
the relativistic approach that these effects are correctly included when the Lorentz invariant 
amplitude (Eq. ([11])) is calculated. 

The difference between the the solid and the dashed curve in Fig. || cannot uniquely 
be separated into relativistic and non-local effects. In order to study the influence of the 
relativistic nuclear structure we performed calculations where we assumed the free relations 
between the upper and lower components of the nuclear wave functions, which led only to a 
slight decrease of the solid curve in Fig. |9|. The difference between the solid and the dashed 
curves in Fig. || is therefore a result of the relativistic, non-local treatment of the production 
operator. In a non-relativistic framework non-local effects have been shown to be important 
in the photoproduction of charged pions in 



B. DWIA 



We now take the final state interaction of the pion with the nucleus into account by 
means of the optical potential described in Sec. |1I E| . The resulting total and differential 



cross sections for 12 C and 40 Ca are shown in Figs. 1C, 11 and 12 in comparison to PWIA 



results. For the case of 40 Ca we took the same parameters for the pion optical potential as 
for 12 C. For E 1 < 230 MeV, i.e. for T^^ 100 MeV, where the absorption is relatively small 
(cf. Fig. D), the distortion of the pion wave function leads to an increase of the total cross 
section, indicating the importance of off-shell effects in the pion distorted wave. For higher 
photon energies the cross section is strongly reduced because of the large imaginary part of 
the pion optical potential. This reduction is stronger in the case of 40 Ca, since the absorption 
of pions increases with the nuclear mass. Besides that, the curves for the two different nuclei 
are very similar in shape, and differ mainly by a global factor. Making an ansatz a ~ A a 
we find a ~ 0.7-0.8 for the case of PWIA. This rather weak A-dependence results from the 
different momentum dependence of the nuclear form factor for 12 C and 40 Ca (cf. Figs. |3] and 
|J): For a given momentum transfer, the 40 Ca form factor is much smaller than the one for 
12 C, which compensates for the fact that there are A amplitudes to be summed for a given 



nucleus (Eq. (|TT|)). For the DWIA results we find a ~ 0.5-0.7, depending on the photon 
energy, which reflects the pion absorption. 

In Fig. the result of a DWIA calculation for 12 C is compared to the A2 data. The 
position of the maximum in the data is well reproduced, but the height of this maximum is 
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underestimated by more than a factor of two. Before we turn to improvements of the DWIA 
calculation, we want to explore two possible sources of uncertainty. 

First one has to be aware that two different but phase-equivalent optical potentials might 
yield the same elastic scattering cross section, but lead to a different behavior of the pion 
wave function in the nuclear interior [[EJ. Since in Eq. (|11|) the pion wave function within 



the nucleus, rather than its asymptotic behavior, is relevant, two phase equivalent optical 
potentials might in principle lead to different cross sections for the coherent photoproduction 
of pions. 

In order to explore this ambiguity, we have also performed calculations using the optical 
potential of |B]] and of |j4"0| , f4"T|| . As can be seen from Fig. |5], these potentials yield pion 
elastic scattering cross sections that are, in the respective energy ranges, comparable to the 
ones obtained with our optical potential. In [13]] the results of a microscopic A-hole model 
calculation for the A-resonance in nuclear matter have been used to construct a position- 
space optical potential for the pion. In Fig. 14" we show the result of a calculation employing 
this optical potential as well as our result in DWIA and PWIA. Both potentials lead to 
a rather similar behavior of the cross section, especially at higher photon energies. For 
the total cross sections, the agreement between calculations employing these two different 
optical potentials is qualitatively the same as for the differential cross section in Fig. WA. 
We have also performed calculations with the parameter sets given in [40] and [IT for T n = 
50 MeV (E~ = 184 MeV). Using these parameter sets, in which the parameters B and C a 
are non-zero, we get a total cross section that differs from the one obtained with our pion 
optical potential by less than 10 % at this specific pion energy. We thus conclude, that our 
empirical parameterization for the pion optical potential accounts for the essential features 
of the pion nucleus interaction in the case of 12 C and 40 Ca. 

As mentioned above, the uNN vector coupling constant and the form factor (Eq. (|5[)), 
which lead to the u contribution in Fig. [7|, are not accurately known. In the elementary 
photoproduction of neutral pions the contribution of the omega graph to the differential 
cross section is about 10 %, depending on the angle. To show the dependence of our result 
for the coherent photoproduction on the omega parameters, we show in Fig. 15 calculations 
employing different values of g^ NN and A. The values g^ NN = 17 and g^ NN = 8 represent 
the upper and lower bound of possible values for this coupling. The resulting uncertainty 
in our results is much smaller than the discrepancy between our DWIA results and the A2 
data. 



C. Medium effects 

By using an optical potential to determine the distorted wave function of the pion, we 
already include the effect of the interaction of the produced pion with the nucleus. The 
intermediate particles contributing to the production operator, however, also interact with 
the nuclear medium. 

For the nucleon propagator we include the resulting medium effects by replacing the 
free propagator with the one obtained from the mean-field in which the bound state wave 
functions were calculated. We thus employ a dressed nucleon propagator Gn that fulfills 

(j/> - m - K7o - V s )G N (p ;x,y) = i5 3 (x-y) . (25) 
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Since Gat is calculated in the mean-field approximation, it contains the static properties 
of a struck nucleon interacting with the residual nucleus. Although it can certainly not 
account for the full dynamics of the struck nucleon interacting with the residual nucleus, it 
is the natural extension of the elementary model to reactions on the nucleus, in which the 



bound and the intermediate nucleon are treated consistently. In Fig. [16] we show results of 
calculations for 12 C using the free as well as the dressed nucleon propagator from Eq. ([25j). 
The increase of the cross section results from the fact, that the total contribution of the 
nucleon terms is the result of an approximate cancellation of direct and exchange graph. In 
the case of the dressed nucleon propagator, this cancellation is more complete, leading to a 
smaller net contribution. Because of the destructive interference of the Born terms with the 
other contributions, this leads to a larger cross section. 

For the A-resonance we can use the A-hole model as a guideline. In flT0| , |TT| it is shown, 
that the medium modifications of the A, in addition to a distortion of the pion, have a clear 
effect on the coherent cross section. It is well known, that the A feels an attractive potential 
of about -30 MeV at normal nuclear density [19|], which has an effect on any production 
process involving the A. Indeed, in ( 3 He, t) and (d, 2 He) reactions a shift of the A-peak 
relative to the reaction on the proton has been observed |]47|j48|| . A considerable part of this 
shift is, however, due to trivial effects e.g. the energy dependence of the nuclear form factor, 
so that it is only partially a consequence of in-medium modifications of the A [49,50]. In 
contrast to that, the A-peak remains essentially unshifted in the photoabsorption on nuclei 
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Calculating the full, relativistic A-propagator in the presence of nuclear potentials is 
out of the scope of this work. We take the attractive self-energy of a A within a nucleus 
schematically into account by shifting its mass in the production operator. In Fig. |16| we 
also show the result of a calculation with the A- mass reduced by 30 MeV (full line). A 
mass shift of the A increases the cross section at lower energies and decreases it at higher 
energies. The value SrriA = —30 MeV has been chosen such that the absolute height of the 
maximum of the A2 data is reproduced (Fig. |1^). Since this number can only be interpreted 
as an average over the entire nucleus, it seems somewhat large in comparison to the depth of 
the A-potential of 30 MeV flTJfl . It must at this point also be kept in mind, that the actual 
value of 5m& needed to fit the data depends to some extent on the value used for of the 
uNN coupling constant. 

Besides the mass, the width of the A is also modified in the medium. A simple analysis 
taking into account Pauli blocking and collisional broadening leads to a width that is rather 
independent of the nuclear density |5^| , but delta- hole calculations yield a A self-energy, 
which corresponds to an increase of the width by as much as 60 MeV [ffl. Including this 
effectively by adding a constant shift of 30 MeV to the energy dependent A-width in Eq. (|J) 
changes the result depicted in Fig. [17] only at energies beyond the maximum, where the cross 
section is rather small. Thus for the A2 data we can clearly distinguish between the effect 
of a modified mass and a modified width of the A. 

This procedure, to take into account the distortions of the produced pion via an optical 
potential and then include medium modifications of the intermediate A by changing its mass 
and its width independently, mimics the essentials of the A-hole model. The qualitative 



agreement with the A-hole model can also be seen in [1G], where the amplitude for the 



coherent photoproduction of pions in the A-hole model is rewritten in such a way, that it 
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contains a distorted pion wave function and a modified A-propagator. In this study it is 
found, that the modifications of the A-propagator, in addition to the pion distortions, have 
considerable effects. 



As can be seen in Fig. [17], the inclusion of a dressed nucleon propagator and a lowered 
A-mass leads to agreement of our calculation with the data around the maximum. At higher 
energies, the situation is not clear. The in-medium production operator seems to lead to 
an underestimation of the data for E 1 > 350 MeV. More experimental data, especially at 
smaller angles, are needed in this energy range before definite conclusions can be drawn. 

Since there is no reliable information available on the properties of an cu-meson at space- 
like momenta in nuclear matter, we did not consider in- medium modifications of the a;-graph. 

The special energy dependence of the differential cross section at a constant angle is 
the result of an interplay between the energy dependence of the production operator and 
the nuclear form factor. In order to circumvent the energy dependence of the form factor, 
and to eliminate the trivial sin 2 6 term (Eq. ([])), we show in Fig. [H| the differential cross 
section for a constant momentum transfer |<f| = 0.1 GeV, divided by sin 2 9 as a function of 
the photon energy for 12 C and 40 Ca. When plotting the cross section this way, the effects of 
in-medium modifications of the production operator are magnified and become visible over 
a large range of energies. Especially an increased A- width now leads to a pronounced effect. 
Thus the new TAPS data, that will cover a large range of angles for each energy [|53|], will 
provide an opportunity to test our method of treating medium effects further. 

In order to also compare to calculations using the A-hole model, we show in Fig. [H] the 
total cross section for 12 C using a medium modified production operator taking all graphs 
into account, as well as that obtained by only including the direct A-term. The result for 
the direct A term agrees qualitatively with results of A-hole calculations. As compared to 
the most recent publications using the A-hole model [§,§J, the height of the maximum in 
the cross section is in our model about the same as in ||, but somewhat above the result 
of |H. The position of the maximum is, however, around E^=220 MeV in [^|,^[, compared 
to 260 MeV in our model. In [|K],|nj the A-hole model is used in combination with an 
empirical multipole decomposition of the production operator, thus including non-resonant 
contributions. Our total result in Fig. [19| agrees with the result of these studies with respect 
to the position of the maximum, but is about 15% higher. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a relativistic and non-local model for the coherent photoproduction of 
pions on spin-zero nuclei, which is applicable from threshold up to the A-region. Since the 
calculation for this process is gauge invariant, and since only a minimum amount of theoret- 
ical input is needed, this reaction represents an attractive field to test our understandings 
of hadrons and nuclei. 

The production operator was derived from an effective Lagrangian for which the free 
parameters were fixed by comparison to the data for elementary processes. Applying this 
production operator to the photoproduction on a nucleus involves an off-shell extrapolation 
with respect to the nucleon. In our model we used the off-shell behavior following naturally 
from an effective field theory. The validity of this extrapolation, can only be tested in 
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comparison to experimental data. Employing a mean field model for the nuclear wave 
function we neglected higher order effects from many-body correlations. 

Our results depend somewhat on the value used for the uNN coupling constant, which is 
uncertain by about a factor of two. Employing the same production operator as in vacuum 
leads to a good reproduction of the shape of the A2 data, but underestimates the absolute 
value by about a factor of two. Agreement with these data could be achieved by taking into 
account medium modifications of the production operator via modifications of the nucleon- 
and A-propagators. 

We conclude that the relativistic DWIA, employing a realistic production operator is an 
adequate approach to the coherent photoproduction of pions on nuclei. More experimental 
data on this process will present a further test of our model and will allow more conclusions 
about the role of medium effects within our model. 
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TABLES 



Nucleus 


v v 


r v 


a v 


V s 


r s 


a s 




(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 




385.7 


1.056 


0.427 


-470.4 


1.056 


0.447 


40 Ca 


348.1 


1.149 


0.476 


-424.5 


1.149 


0.506 



TABLE I. Strengths, reduced radii and diffusivities for the relativistic scalar and vector 



mean-field potententials, respectively. 
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FIGURES 




FIG. 1. Feynman diagrams contributing to the photoproduction of neutral pions. 




19 





20 




FIG. 5. Fit results for elastic pion scattering as described in the text. The data are taken 
from Ref. ^2|. The dashed lines show results for the parameter set C from |^(| at 50 MeV and for 
the optical potential from 44] at 157 and at 220 MeV. 
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FIG. 6. Energy dependence of the s- and p-wave parameters for the pionic optical potential 
resulting from our fit as described in the text. 
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FIG. 7. Contributions of the different fundamental diagrams to the total cross section for 
C(7,7r°) 12 C in PWIA. 
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FIG. 9. Total cross section for 12 C(7, 7r°) 12 C in PWIA resulting from the direct A-graph (solid 
line) in comparison to non-relativistic, local calculations as in [ p^Jl~7[| (dotted line) and including 
the kinematical correction from (dashed line) as described in the text. 
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FIG. 11. Differential cross section in DWIA and PWIA for 12 C(7,tt° 
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FIG. 13. Comparison of the A2-data for 12 C to a DWIA calculation. 
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FIG. 14. Comparison of the A2-data for 12 C to DWIA calculations employing two different 
optical potentials and a PWIA calculation. 
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FIG. 15. Comparison of the A2-data for 12 C to DWIA calculations employing different loNN 
couplings and cutoffs. The solid curve was obtained with the values used throughout this work. 
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FIG. 16. Total DWIA cross section for 12 C(7, 7r°) 12 C using the free production operator (full 
line), using the dressed nucleon propagator (dashed line) as described in the text, and with a mass 
shift of the A 5m A = -30 MeV (dotted line). 
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FIG. 17. Comparison of the A2-data for 12 C to a DWIA calculation employing a 
medium-modified production operator with a dressed nucleon propagator and a shifted A-mass 
(Jm A =-30 MeV). 
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FIG. 18. Differential cross section divided by sin 2 9 for 12 C( 7 , tt°) 12 C and 40 Ca(7, 7r°) 40 Ca for a 
fixed momentum transfer of 0.1 GeV: free production operator (dashed line), in- medium production 
operator with a dressed nucleon propagator and a A-mass shift 5mA =-30 MeV (full line), and with 
an in-medium production operator including additionally an increased A-width 5T& = 30 MeV 
(dotted line). 
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FIG. 19. Total cross section for 12 C(7, vr°) 12 C using an in-medium production operator 
(8m a =-30 MeV, 5Fa = 30 MeV), for a complete calculation, as well as for a calculation where 
only the direct A-graph is included. 
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